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Abstract 

Clusters of infected individuals are denned on data from health laboratories, but this 
quantity has not been defined and characterized by epidemy models on statistical 
physics. For a system of mobile agents we simulate a model of infection without 
immunization and show that all the moments of the cluster size distribution at the 
critical rate of infection are characterized by only one exponent, which is the same 
exponent that determines the behavior of the total number of infected agents. No 
giant cluster survives independent on the magnitude of the rate of infection. 
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1 Introduction 



In a small and highly urbanized nation like Singapore dengue outbreaks or epi- 
demics are identified as "clusters" . A dengue cluster or focus of transmission 
is defined as at least two confirmed cases, with no recent travel history, that 
are located within 200 m of each other (taken as the flight range of the Aedes 
aegypti) and whose dates of the onset of symptoms are within three weeks of 
each other [1]. Some efforts have been directed towards the characterization 
of 'SIS' models of infections, or epidemics without immunization [2,3,4], that 
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is the state of the particles are healthy or infected, and are susceptible to re- 
infection after healing, thus the name of the model (SIS: susceptible-infected- 
susceptible). Analytical and numerical expressions describe the dynamics of 
the SIS model in terms of the rate of spreading A, the evolution of the survival 
probability of infection P(t), the mean number of infected agents n(t) and the 
mean square distance of spreading R 2 (t) in time, which are quantities difficult 
to compare with real data of epidemics. This work suggests an application of 
potential comparison with public health data, analyzing a scaling function for 
clusters numbers on a SIS model of infection. 

The second important ingredient of this work is the mobility of agents, con- 
trasted to most of the models of epidemy where the population is modeled by 
static networks [2,3,4]. In a previous work [5] we showed that the SIS model 
of infection on a system of mobile agents has critical exponents which depends 
on the density of the system, i.e spatial correlations and mobility of the agents 
play an important role. We obtained a crossover from mean field behavior for 
low densities to static 2D-lattices for higher densities. Here we use our model 
of mobile agents to define clusters of infections and analyze its dependency 
on the rate of infection A (defined in detail bellow) and on the mobility of the 
agents. 

We propose a time-evolving network model: A link between two moving agents 
is created when they collide with each other and there is transmission of the 
infection among them (i.e through infected-susceptible interactions), the link 
lasts a characteristic time of infection. 

We find that the network of clusters of infections remains disconnected and 
no matter how large the rate of infection, no giant cluster is formed. We show 
that in the transition to spreading, the moments of the cluster size distribution 
are described by an exponent (3, which is the exponent that characterizes the 
fraction of infected mass Fjm = Ninf/N, defined as the ratio of the number 
of infected agents (N In f) and the total amount of population (N). Thus the 
number of clusters depends on A, and mobility and spatial correlation of the 
agents influence its dependency. 



2 Model 



N soft disks, with radius ro = 1, represent agents which move in a two dimen- 
sional cell of linear size L, with density p = N/L 2 . The system has periodic 
boundary conditions and is initialized as follows: the agents are placed in the 
cell with the same velocity modulus v and randomly distributed directions, 
positions and states: 'infected' or 'susceptible'. If a susceptible agent i collides 
with an infected agent j (i.e |t*» — rj\ <= 2vq), then % becomes infected. Each 
infected agent heals and becomes susceptible again after a fixed number of 
time steps, called the 'time of infection' (At in f), which is a free parameter of 
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the model. 

The physical interaction of the agents is modeled by molecular dynamic with a 
leap-frog integration method [7], the interaction potential is a 12 — 6 Lennard- 
Jones truncated potential (see more details in [5]). 

The resulting model is a contact process [6] , where the infected species become 
extinct unless the infection spreads rapidly enough. The transition between 
survival and extinction depends on a critical rate of spreading A c that marks 
the transition into an absorbing state. The infection rate A is defined as the 
number of agents one agent infects before healing. For this system, 

A = At inf /T f , (1) 

where Tf is the characteristic time of flight between two collisions, which is 
determined by the density (p) and the mean velocity of agents ((v)). The 
critical exponents of the transition to spreading were presented by us for the 
same kind of system [5], where the study was done in terms of the fraction of 
infected individuals (Fjm)- Here we go further and characterize the behavior 
of the clusters of infected individuals. When agent j infects agent i a link is 
created among them, the link lasts until one of them heals, meanwhile each 
of them continues making links with other susceptible agents through the 
same rule. A cluster is thus defined as a group of infected agents connected 
by links. Note that in contrast to percolation, where clusters are given by 
occupied lattice sites connected by nearest-neighbor distances, for this model 
each cluster gives a group of agents infected in a given period of time linked 
by a relation of contagion. Isolated infected agents are regarded as clusters of 
size unity and any cluster consisting of s connected agents is an s — cluster. 
We borrow the notation from Stauffer's book on percolation theory [8] and 
define here n s = N s /N as the number of s-clusters per agent, where N s is the 
number of clusters of size s and N the total number of agents in the system. For 
different values of A, in the next section we present the results of the calculation 
of the first three moments of the cluster size distribution. Namely: J2 s n s, 
J2 S sn s, J2 S s 2 n s . Those quantities give us, respectively, information about: the 
total number of clusters, the fraction of infected agents and the mean size of 
clusters. In order to keep the analogy with the calculation on percolation, we 
sum over all values of s excluding the largest cluster (S ma j or ). We also present, 
the calculations of F ma j or = S ma j or /N, the fraction of agents that belong to 
the largest cluster and Fjm = N in f/N, the fraction of agents that are infected. 



3 Results 

For a fixed density, we vary A (Eq. 1) changing the time of infection (At in f). 
Starting with half of the population infected, for rate of infections near A c , a 
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Fig. 1. Left: Fraction of infected individuals from surviving trials versus time at 
A = A c , starting with half of the population infected. At the top, the results for 
p = 0.05 and A = 1.06 and at the bottom p = 0.46 and A = 0.68. System sizes 
N = 32 x 32, 64 x 64, 128 x 128 (from top to bottom). Right: Quasi-stationary 
fraction of infected agents versus A for the same densities (Top: p = 0.05. Bottom: 
p = 0.46). 

given trial may end into the absorbing state after a few time steps or it may 
survive fluctuating with a quasi-stationary fraction of infected agents, marked 
with windows in the left-side of Fig 1. The calculations are made averaging 
on time at the quasi-stationary state, which is described by the surviving 
trials following an initial transient. The number of time steps of this transient 
depends on A and on the system size L (see left side of Fig. 1). The data here 
illustrate how the mean fraction of infected agents Ff^it) (the superscript 
denotes an average restricted to surviving trials) approaches its stationary 
value Fjm{^,N) (in the following, we write Fim{^,N) just like -Fjm(A)). In 
the right side of the same figure we see the graph of Fjm(A), which becomes 
sharper increasing the system size. We analyze in detail the number of clusters 
for the two density values p = 0.05 and p = 0.46, which have critical rate of 
spreading A c = 1.06 and A c = 0.68 respectively. Note that at the critical 
density A c , surviving trials tend to stationary values only in the limit L — > oo. 
The left side of Fig. 2 is only for pedagogical reasons, in order to illustrate how 
the number of clusters looks in the quasi- stationary state, we see snapshots of 
the clusters of infections for different systems densities and the same rate of 
infection A = 1.5, here N = 10 x 10. 

For p = 0.05 and A G [1, 200], the right side of Fig. 2 shows the variation of 
-P/m(A) and N = 32 x 32 averaged over 20 different realizations. The insets 
show the change in time of F ma j or and J2 S n s, for only one realization with A = 
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Fig. 2. Left: Snapshots of cluster sizes of infected agents for systems with different 
densities: (a)p = 0.05, (b)p = 0.20, (c)p = 0.40 and (d)p = 0.80, in all cases A = 1.5. 
Right: Quasi-stationary fraction of infected agents varying A over three orders of 
magnitude (Average over 20 realizations for p = 0.05 and N = 32 x 32). The insets 
show the fraction of infected agents in the largest cluster (lower value) and the the 
first moment of the cluster size distribution (upper value) vs. time, at A = 1.08, 
A = 10.0 and A = 108.0. 
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Fig. 3. First three moments of the cluster size distribution , fraction of agents in the 
largest cluster (F ma j or ) and fraction of infected agents {Fjm) vs. A. Average over 
50 trials, system size N = 64 x 64. Left:p = 0.05. Right:/? = 0.46 

1.08, A = 10.0 and A = 108. In contrast to percolation results, in this model 
there is no significant variation of F ma j or with A, and the relation F ma j or <C 
Fjm remains. Moreover, the number of clusters J2 S n s grows considerably only 
near A c . 

In Fig. 3 for p = 0.05 and p = 0.46, we plot the behavior of the cluster numbers 
near their respective A c . As the largest cluster remains small compared to the 
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Fig. 4. Same results of Fig. 3 plotted vs. (A — A c ). The solid lines are regressions of 
the form mj(A — Ac)' 3 with rrtj the coefficient of the ith moment. Left: A c = 1.06, 
(3 = 0.66, mo = 0.321, m\ = 2mo, and ni2 = Qm,Q. Right: A c = 0.68, (5 = 0.56, 
niQ = 0.386, mi = 2.3mo, and mi = 7.5mo 

total number of agents (S ma jor *C N), we have i 7 Vjw(A) ~ J2 S sn s- Additionally 
one can see that J2 S sn s and J2 S s<2n s show the same critical behavior as -Fjm (A), 
plotted in detail in Fig. 4. We observe that all the moments of the cluster 
size distribution present exactly the same critical behavior than F IM , namely 
(A — A,.)' 3 , where (3 depends on the density of the system. 



4 Conclusions 



This work showed that the cluster size distribution of infected individuals 
is described in terms of the spreading rate (A) and the same exponents (j3) 
previously known for the total mass of infection. Although the agents are free 
to move there is a homogeneous size distribution of infected clusters at the 
critical rate of infection, and we did not find any critical exponent associated 
with the cluster sizes. Comparing with the traditional SIS model on a static 
network we confirm that mobility and spatial correlations change the value of 
the critical exponent f3 of the fraction of infected population, and to the same 
extent the cluster size distribution of infection. 
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